Mesoscopic anisotropic magnetoconductance fluctuations in ferromagnets 
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The conductance of a ferromagnetic particle depends on the relative orientation of the magne- 
tization with respect to the direction of current flow. This phenomenon is known as "anisotropic 
magnetoresistance" . Quantum interference leads to an additional, random dependence of the con- 
ductance on the magnetization direction. These "anisotropic magnetoresistance fluctuations" are 
caused by spin-orbit scattering, which couples the electron motion to the exchange field in the fer- 
romagnet. We report a calculation of the dependence of the conductance autocorrelation function 
on the rotation angle of the magnetization direction. 

PACS numbers: 75.75.+a, 72.15.-v, 72.25.-b, 73.22.-f, 75.30.Gw 
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One hallmark of phase-coherent transport is the phe- 
nomenon of "universal conductance fluctuations" , ran- 
dom, but reproducible variations in a sample's conduc- 
tance as a function of the applied magnetic field or the 
Fermi energy P,Q>0,0,01- The magnitude of the conduc- 
tance fluctuations is of order unity, in units of the con- 
ductance quantum e 2 //i, and does not depend on specific 
sample properties, such as the impurity concentration, 
the meterial, shape, or method of preparation. 

Recently there has been both theoretical and experi- 
mental interest in mesoscopic transport in itinerant ferro- 
magnets. The experimental interest stems from the abil- 
ity to fabricate ferromagnetic conductors small enough 
that transport through the magnet is predominantly co- 
herent |a, [3. The theoretical interest is motivated by 
the rich variety of ways through which random impurity 
scattering can affect the properties of an itinerant fer- 
romagnet. Theoretical predictions exist for the effect of 
domain walls on weak localization and conductance fluc- 
tuations |a |9( as well as for the combined effect of spin- 
orbit interaction and impurity scatte ring on weak local- 
ization and magnetic anisotropy [flj. Although dis- 
ordered ferromagnetic conductors display different phe- 
nomena than their normal-metal counterparts, the theo- 
retical framework to describe them is rather similar. In- 
deed, the methods of diagrammatic perturbation theory 
developed for electron transport in disordered metals can 
be applied to ferromagnets by modifying the single par- 
ticle Hamiltonian taking into account the exchange field 
and/or spin-orbit interactions. 

In this Letter, we address the mesoscopic contribu- 
tion to a ferromagnet's anisotropic magnetoresistance. 
Anisotropic magnetoresistance is the phenomenon that 
a magnet's resistance depends on the orientation of the 
magnetization resulting from a combination of spin-orbit 
coupling and orbital magnetic effects For a single 

domain magnet, the resistance is a smooth function of 
the magnetization direction. The mesoscopic effect de- 
scribed here consists of an additional and faster random 
dependence on the magnetization direction that is differ- 
ent for each sample, but reproducible for a given sam- 



ple. This situation is not very different from the case 
of standard universal conductance fluctuations in a nor- 
mal metal, where the random magnetic-field dependent 
fluctuations are superimposed on a systematic magneto- 
conductance. 

There are two possible mechanisms through which the 
magnetization direction can affect the interference cor- 
rection to the conductance. First, a change of the mag- 
netization direction causes a change of the internal mag- 
netic field, which directly affects the orbital motion of the 
electrons via a change of Aharonov-Bohm phases. Sec- 
ond, a change of the magnetization direction causes a 
change of the exchange field, which affects the motion 
of the electrons via spin-orbit scattering. The first ef- 
fect would be dominant if the magnetic flux through the 
cross-section of a phase coherent volume is of the order 
of the flux quantum. For many magnetic materials, the 
phase coherent lengths can be small and this effect can 
be neglected (see discussion in Ref . ITo(l . In what follows, 
we assume that this condition holds, and that the second 
effect dominates the mesoscopic anisotropic magnetore- 
sistance. For the same reason, we ignore any effect of an 
applied magnetic field used to change the magnetization 
direction. 

We consider an ensemble of ferromagnetic particles, 
each with a different configuration of impurities and cal- 
culate the conductance autocorrelation function 



C(9) = (G(m)G(m')) - <G(rh)) : 



(1) 



where is the angle between the magnetization directions 
rh and rh' and the brackets (. . .} denote the ensemble 
average. The vectors rh and rh' are defined to have unit 
length. The Hamiltonian for a ferromagnet with spin- 
orbit scattering is 
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2m 



a/3 



V a /3 



(2) 



where a and /3 are spin indices, a z the Pauli matrix, 
the magnetization direction rh is taken as the spin quan- 
tization axis, and Ez = Mb-B cx is the Zeeman energy 



2 



corresponding to the exchange field B ex . We perform 
the ensemble average at a fixed chemical potential \x and 
exchange field i? GX , rather than at self-consistently de- 
termined /i and B cx . Although the omission of the self- 
consistency conditions is known to affect average d q uan- 
tities, it is believed not to affect fluctuations [la. ll4L fl5| . 

The random potential V in Eq. describes the effect 
of elastic impurity scattering and spin-orbit scattering, 
respectively. Its Fourier transform is 



Vak,,3k' — Vk-k' 

-C k -((k'xk).K 



(3) 



ei(T 



e 2 0- y )a/3, 



where ei and e 2 are unit vectors perpendicular to each 
other and to m such that ei x e 2 = m. The random 
potentials V and V so are assumed to be uncorrelated 
and Gaussian white noise, with r.m.s. strength v and v so , 
respectively, 

(VM = V 2 S(q q'), (V*°V°?) = v 2 s J( q - q'). (4) 

In the leading order Born approximation, the scattering 
time T a for spin-independent impurity scattering of elec- 
trons with spin a is given by 
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2-KVlTl 



(5) 



where v a is the density of states of electrons with spin a. 
Similarly, for spin-conserving and spin-flip scattering off 
V s °, one has the mean free times 
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respectively, where kp a is the Fermi wavevector for spin 
a electrons. In a realistic ferromagnet, the kinetic energy 
and the random potential will not have the simple form 
assumed in our calculation, which implies that the rela- 
tionships between the scattering times implied by Eqs. 
(J5J and 10 need not hold. Although we use the simple 
model described above to set up our calculation and to 
define the scattering times, these are then considered in- 
dependent for the rest of the calculation [except for the 
equality in the second line of Eq. © , which follows from 
detailed balance]. 

Throughout the calculation, we assume that r <C tii , 
tj_. This implies that all Green functions appearing in 
intermediate phases of the calculation can be averaged 
over all directions of the momentum. We also assume 
that phase coherence is preserved over the entire sample. 
In a sample with size L larger than the phase coherence 
length L^, our answer would be modified as C(8,L) ~ 
C{6,L ( f ) ){L ( f ) /L'). In this case, the angle over which the 
conductance typically fluctuates is then determined by 
Lj, instead of L. 
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FIG. 1: Dyson Equation for Diffuson ladder. The dotted line 
indicates a scattering event. 



We now describe the details of our calculation. For the 
retarded Green function C? R , averaged over the random 
potential and over all directions of the momentum, we 
find 

k, m)}- 1 = u - e a (k) + -L + -L- + (7) 

Z Try L Try 1 1 Try I 



where s a (k) = hrk /2m — fj, 



is the energy of 



an electron with spin a and momentum Kk. In order to 
calculate the conductance autocorrelation function |JT]|. 
we need to consider the Diffuson and Cooperon propa- 
gators of diagrammatic perturbation theory. Again, in 
view of the inequality r <g. T|| , t±, we only need Diffuson 
and Cooperon propagators averaged over all momentum 
directions. Since the Cooperon and Diffuson propagators 
are related by time reversal, 



(8) 



it will be sufficient to calculate the Diffuson only. 

The Diffuson propagator is defined by the ladder dia- 
grams shown in Fig.^ The solid arrows in Fig^denote 
the impurity-averaged Green functions J7J) . The two legs 
of the ladder refer to the two magnetization directions 
m and m'. For both magnetization directions we use 
the convention that the magnetization direction is the 
spin quantization axis. This is the natural choice for 
ferromagnets: Since E%t > 1 in a typical ferromagnet, 
with this convention only ladder diagrams for which the 
spin indices of retarded and advanced Green functions 
are pairwise equal at all times need to be considered; 
contributions with different spin index for retarded and 
advanced Green functions dephase within a mean free 
time and do not contribute to the Diffuson propagator. 
One should note, however, that this convention implies 
that the directions of "spin up" and "spin down" in the 
upper and lower legs of the ladder correspond to different 
physical directions if m ^ rh'. 

Summing the ladder diagrams of Fig. ^ we then find 
that the Diffuson obeys the 2x2 matrix equation 

1 
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K ai D(oj, q, 9) lf 3 



a f3 



2irv a T a 



(9) 



Here K is a 2 x 2 matrix, with diagonal elements given 
by 

2 1 - cosfl" 



(10) 
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FIG. 2: Leading diagrams for conductance correlator. The 
wavy lines represent the current vertex J and the shaded box 
represents either the Cooperon or Diffuson propagator. 

where D a = Wp Q T a /3 is the diffusion constant. The off- 
diagonal matrix elements contain a phase factor that de- 
pends on the precise choice of coordinate axes perpen- 
dicular to m and rh', cf. Eq. J2J|. In all final expressions, 
the off-diagonal elements of K only enter through their 
product, which is independent of this choice, 

K n K u = _3IL(1 + cos Of. (11) 

Once the Diffuson is known, the Cooperon is calculated 
via Eq. (jSJ. For the special case 8 = 0, the result for C 



was previously obtained by Dugaev et al. 10]. 

We can now proceed to calculate the conductance cor- 
relation function C{9). We are interested in the conduc- 
tance correlations at zero temperature, which allows us 
to set u> = in our expressions for the Diffuson and 
Cooperon propagators. We consider a coherent rectan- 
gular sample with sides L x , L y and L, , with a current in 
the z direction. Following Refs. |2> y, we then find that 
the conductance autocorrelation function is given by 

C(0) = ]T tr [JD(q) JD(q) + JC(q)JC(q)], 

z q 

(12) 

where the current vertex reads 

47T 

Ja(i = y VFa T a<W (13) 

and the vector q is summed over the values q x = irn x /L x , 
q y = Tmy/Ly, and q z = nn z /L z with n x , n y = 0, 1, 2, . . . 
and n z = 1,2, . . .. Without the prefactor 3/2, Eq. JHJ) 
is precisely the contribution from the diagram shown in 
Fig. El The factor 3/2 in front of Eq. (jT2"|> accounts for 
other diagrams that contribute to the conductance fluc- 
tuations, whose net contribution is 1/2 times that of the 
diagram of Fig. H 0,0]. Substituting our results for the 
Diffuson and Cooperon propagators, we find 



q ± 



((L z q/n) 2 + a ± (6)) 2 {{L z q/ir) 2 + a±(7r - 9)) 2 



(14) 



where 



a±(0) 



i , i , T n E i+ T i\\ E i 



(1 - cos 0)) 
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/ (1 + cos6») 2 
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(1 - cos 6) 



(15) 



and E a = D a (-K/L z ) 2 is the Thouless energy for spin 
a. Note that the parameter that governs the importance 
of spin-orbit scattering is the product r a j_E a or T a \\E a , 
which is the ratio of the spin-orbit time and the Thouless 
time, which is the time to diffuse through the sample. 

The expression for a±(8) simplifies in two limiting 
cases. If 6 = 0, one has a + = 2/t^±E^ + 2/t^£^ 
and a_ = 0, showing the presence of universal conduc- 
tance fluctuations in a ferromagnet. The correspond- 
ing eigenvalues for the Cooperon contribution are found 
by setting 9 = tt, a+(7r) = 2/t^±E^ + 2/t^\E^ and 
a~(ir) = 2/ti±Ei + 2/rniE^. Another simple limit is 
that of a half-metal, a ferromagnet with vanishing den- 



sity of states for the minority spins. For a half metal, 
the only relevant time and energy scales are the scat- 
tering time Tf|| for spin-preserving spin-orbit scattering 
of majority electrons and the majority electron Thouless 
energy E^ . One then finds that only one root a± is rele- 
vant, a{8) = (1 — cos 6>)/-r-f || E^. 

The sum over wavevectors in Eq. 1141 can be performed 
analytically for a quasi one-dimensional sample. Setting 
n , 



n y = in the summation, one finds 



C ( e ) = E [^V^W) + F(nVa ± (7r - 0))] ,(16) 
where F(x) = 3e 4 (— 2 + ircothx + x 2 sinh -2 x)/2x i h 2 . 
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FIG. 3: The correlation function of the conductance at dif- 
ferent directions of the magnetization, for various strengths 
of the spin orbit scattering. Results shown here are for a half 
metal with cubic geometry. 



Note that for 9 = 0, Eq. (|16f) reproduces the known re- 
sults varG = (e 2 //i) 2 (l/15) for strong spin-orbit scat- 
tering and varG = (e 2 //i) 2 (4/15) for weak spin orbit 
scattering. For quasi 2D and 3D samples C{9) can be 
computed numerically. The dependence on the spin- 
orbit scattering is qualitatively similar for all these cases. 
Shown in Fig. [31 is C{9) for a half metal with L x = 
L y = L z . The top dashed line in Fig. |3]is the variance 
of the conductance in the absence of spin-orbit scatter- 
ing. Without spin-orbit scattering, there is no angle- 
dependent mesoscopic correction to the conductance, so 
C{9) is independent of 9. For t^\\ <C 1/Ef, conductance 
fluctuations saturate at half their value without spin- 
orbit scattering. Changing the magnetization by a small 
angle 9 C changes the mesoscopic conductance correction 
enough to lose all conductance correlations. Our calcu- 
lation shows 



(ni^t) 1 / 2 ~ (r n M)^l/L 



1/2; 



(17) 



where / is the mean free path. In a realistic ferromagnet, 
the quantitative form of C(9) is different, although the 
qualitative picture, including the estimate for the cor- 
relation angle 9 C is the same as for the half metal (see 
Eq.HSJ. 

Let us estimate the correlation angle 9 C for the spin- 
orbit induced mesoscopic conductance fluctuations. For 
the highly disordered ferromagnetic wires used in the ex- 
periments of Refs.0,0, the mean free path I is of the order 
of a few nm. Taking the spin-orbit times tm and tj_ within 
an order of magnitude of the elastic scattering time r (as 
is appropriate for Co 0), we find 9 C ~ (1 x lCr 8 m)/L. 
(Recall that L has to be replaced by the phase coher- 
ence length Lcj, if < L.) This would be sufficiently 
small to explain the few conductance oscillations seen in 



the experiment of Ref. 0, for which ~ 30nm and the 
conductance was measured as a function of an external 
magnetic field that changed the magnetization direction. 

It is instructive to compare the correlation angle 9 C for 
spin-orbit induced conductance fluctuations considered 
here to the correlation angle arising from the coupling of 
the electron's charge to the internal magnetic field. The 
latter is ~ &o/&, where <£> is the magnetic flux through 
the sample and $o is the flux quantum. Taking the in- 
ternal magnetic field to be ~ 2T, as is appropriate for 
Co, one finds a correlation angle ~ (2 x 10 _15 m 2 )/L 2 . 
Hence, with the parameters taken above, the orbital ef- 
fect will dominate for samples with size L > 2 x 10 _7 m. 
This is in agreement with Ref. 7, where it was shown that 
the orbital effect alone cannot account for the observed 
conductance fluctuations 0. 
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